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zation problem. Although many diverse effects are not considered, 
the results nevertheless provide a guide to limits of performance of 
plane fin/heat-pipe radiators using recently developed materials. 
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INTRODUCTION 

Because of its importance for space power and thermal manage- 
ment, heat rejection in space radiators has received considerable 
attention over the years. Basic results of heat rejection studies have 
been reported in texts on radiation (Sparrow and Cess, 1978; 
Edwards, 1981). Minimum mass designs were addressed in 1960 
(Bartas and Sellers, 1960); more recently, comparisons have been 
made between pumped-loop and heat-pipe radiators (Furukawa, 
1981), with finned heat-pipes being the choice today. 

Despite the extensive previous work, much of which is summa- 
rized by Furukawa (1981), there is still interest in the problem as 
new materials are introduced (McDanels et a!., 1991) and as larger 
power plants are considered for space applications. In particular, 
new anisotropic composites, where the transverse thermal conductiv- 
ity may be an order of magnitude less than along the fiber axis 
(Rhee et ah, 1991; Rovang et al., 1991; Materials and Science Cor- 
poration, 1986), have raised questions about anisotropic conduction 
effects in the radiation fins. Furthermore, despite the extensive use 
of heat-pipe radiator designs, the effective heat- pipe/ radiator- fin 
interaction has not been addressed, with previous analyses using the 
standard specified root-temperature boundary condition. 

The present paper seeks to address these shortcomings, as well 
as to investigate several minimum-mass heat-pipe radiator designs. 
Two-dimensional, anisotropic heat conduction equations are formu- 
lated for plane fins with radiative conditions on the surfaces. These 
are reduced to axial integration using integral techniques (Kakac 
and Yener, 1985), with anisotropic effects contained in a surface 
radiation Biot number. Effective fin/heat-pipe boundary conditions, 
and a shape factor are introduced following previous examples 
(Lund, 1986, 1989). The axial integration was carried out numeri- 
cally, with approximately 100 steps along the length of the fin, and 
fin performance curves were derived. Typical designs were analyzed 
for minimum mass using the performance curves, and an optimum 
fin length, for various fin efficiencies and for heat-pipe vapor 
temperatures of 400 and 800 K. It was found that an efficiency of 
0.5 yields the minimum mass, which agreed with previous results 
(Sparrow and Cess, 1978; Bartas and Sellers, 1960); however, it is a 
very fiat minimum, and other efficiency values can be used at the 
designer’s option without incurring significant penalty. 


The thermal analysis of the space radiator considered here may 
be divided into the following parts: (a) radiation-enclosure analysis 
of the exterior to the radiator fin, or panel; (b) conduction analysis 
internal to the fin; (c) analysis of the heat-pipe condenser adjoining 
the panels; and (d) the interaction between these effects. Here the 
objective is to separate each effect, as much as possible, by the in- 
troduction of fin efficiency, effective surface conductances, and shape 
factors. This simplifies the model equations, as well as generalizes 
the results to a wider class of designs, as shown previously (Lund, 
1986, 1989). 


Exterior Radiation Model 

We consider an element of radiator fin, as shown in Fig. 1, 
where the net heat rates (per unit width) leaving the element from 
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Figure 1 . — Control volume for integral analysis of fin. 


surfaces 1 and 2 are dqj , and dq 2 , and where the corresponding 
surface temperatures at the location X are T bl and T b2 . Assum- 
ing these surfaces to be gray (diffuse with constant emissivities), 
and exchanging radiant energy with surrounding gray surfaces (j at 
temperature T b j), the net heat rate from one side of the element 
(say, side 1), may be stated as 

d <u - H M T bi' T i) dA (i; 

where ^"is the usual gray-body view factor, or transfer function 
(e.g., Siegel and Howell, 1981; Edwards, 1981), which satisfies the 
summation 


Even at the higher temperature, the Biot number was found to 
be a very small number, such that anisotropic effects are negligible 
for the caseB considered. Overall, conduction affects the minimum- 
mass designs as a weak function, as compared to the material den- 
sity, and carbon/carbon composites had the lowest mass of the 
designs considered. The effect of the interface shape factor was also 
evaluated. 


£j *1} m ‘l (2) 

Combining Eqs. (1) and (2) results in the exterior radiation model 

d <h' "‘.'(Tbl - T rl) dA (3) 


Minimum-mass radiator designs have particular importance for 
overall space power system design, as was indicated for solar dy- 
namic power (Lund, 1990). However, radiators in actual systems 
operate under many constraints not considered here; in this sense, 
the present analysis may be considered as an unconstrained optimi- 


where the surface- 1 effective radiation temperature of the surround- 
ings is 
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with similar results for surface 2. Thus, the external radiation 
problem reduces to determination of the T r ’s (which may have 
some X-dependence through the transfer functions). 

In the space environment, however, the effect of the radiation 
temperature is small. Even if the heat-pipe condenser temperature 
is as low as 600 K and the surrounding radiation temperature is as 
high as 300 K, this is a ratio of only 0.5 for which there is negligible 
effect on the overall rate of heat transfer, or fin efficiency (Sparrow 
and Cess, 1978; Edwards, 1981). Therefore, for the remainder of the 
present analysis, T r will be neglected, which still accounts for most 
space radiator problems of interest. For ratios greater than 0.5, 
corrections to the results of the present work can be made with 
reference to the efficiency curves of Sparrow and Cess (1978). The 
effect of neglecting T r is that the fin conduction problem becomes 
symmetric, regardless of asymmetries in the radiation environs. 

Conduction Fin Model 

With reference to Fig. 1, an energy balance on the control vol- 
ume yields 

‘lx “ dc ll + qx+dx 


or, with Eq. (3), 
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where the rate of heat conducted is the integral of the flux over 
the half cross section. With the axial heat flux related to the axial 
temperature gradient by the axial conductivity k x , Eqs. (5) and (6) 
combine 
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On the other hand, the surface heat flux depends on the transverse 
conductivity k y , such that the temperature in Eq. (7) must satisfy 
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and the symmetry condition 
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Y=0 


Based on the assumption that there is negligible heat loss from the 
tip of the fin, the axial boundary conditions are the insulated tip 
condition 


dT 

ax 


( 10 ) 


X- 0 


and an equivalent convection-type condition at the heat-pipe 
interface 
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where h e and T a are effective parameters (to be determined) 
associated with particular heat-pipe designs. 


Now, introducing the scaled coordinates and temperature 
x = X/L, y = Y/b, 0 — T/T e , Eqs. (7) and (8) may be stated as 
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where the average and surface temperatures are 
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and 
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and where the radiation/conductance and Biot numbers are 


N - <r<rT*L 2 /b f k x 


(12) 

(13) 

(14) 

(15) 

( 10 ) 


and 


P - e*T*b f /k y (17) 

The effects of transverse and anisotropic conduction are 
contained primarily in the surface Biot number parameter. These 
effects have more significance and the parameter becomes larger at 
higher temperatures and larger fin thicknesses. 

A temperature profile which satisfies the centerline symmetry 
condition and Eq. (15) is given by 

tf(x,y) = 0 c (x) - K(x) - tf,(x)]y 2 (18) 

where $ c is the centerline temperature. Substitition of Eq. (18) into 
Eq. (13) yields 

*c “ 0 . * (19) 

and substitution of Eqs. (18) and (19) into Eq. (14) yields the 
average temperature 


( 20 ) 


Differentiating Eq. (20) twice and substituting into Eq. (12) then 
yields the final equation for the surface temperature 
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It is seen that 5 is always higher than by some amount, but 
that for vanishingly small surface Biot numbers, & « B % , and 
Eq. (12) reduces to the conventional radiation fin equation. With 
the integral approximation, Eqs. (7) or (12), the boundary condition 
(Eq. (11)) cannot be satisfied for all y. 


Here we take it to hold in the mean 
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where the effective end-condition number is 
1 - h„L/k x 
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(23) 



Results of Thermal Model 


An example of the numerical integration of Eq. (21), with 
Eqs. (20) and (22), is shown in Figure 2. The end-condition 
parameter 7 (the heat-pipe interface condition) has considerable 
effect on the solution; the limit 7 — ► «o corresponds to the usual, 
specified root-temperature condition. The difference in surface and 
average temperatures is also noted for a surface Biot number of 0.2, 
which accounts for anisotropic conduction effects. 



Figure 2.— Typical average and surface temperature 
distributions in radiator fin for N = 1 and p * 0.2. 


Of particular interest is the overall rate of heat transfer through 
the fin, 

L 1 

q f = J q dX * L eo j T*dx 

0 0 

or 

q f = L caT* (24) 

where the fin efficiency is give by 

1 

- jVdx ( 25 ) 

0 

This function is shown in Fig. 3. For f) = 0.2, which here 
accounts for anisotropic conduction, there can be a dramatic effect 
on efficiency for small N. 

From a design and optimization point of view, the radiation/ 
conductance number N is really an unknown since it includes the 
(as yet) undetermined length or thickness of the fin. Therefore, it is 
advantageous to invert Eq. (25) and express N as the function 
N = N(7,i7 f ,/9). This was accomplished numerically in a double 
iteration scheme using the simplex algorithm, with results as shown 
in Fig. 4. 

Essentially, for a given fin length and efficiency (and Biot 
number), N is the inverse fin thickness, (seen the definition, 

Eq. (16)): 

b f > e,T®L 2 /k x N(h e L/k x) (26) 

Thus, for a selected fin efficiency and Biot number, a curve-fit 
of Figure 4 may be used to determine bf directly (and iteratively if 
ft is not known, but is to be calculated). 

Heat-Pipe Condenser Model 

A diagram of a typical heat-pipe condenser section is shown in 
Fig. 5, where the thicknesses of the pipe and liquid wick are b p and 



Figure 3. — Radiator fin efficiency function. 
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Figure 4. — Radiator fin design chart. 
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( 31 ) 



* “ e p < 7 T*(b p /k p ♦ b,/k,) 

For the range, 0 < 6 < 1, Eq. (30) is well represented by 

, p = 1 - 0.736i(°- 252 -° 0868lni ) (32) 

With these definitions, the total heat transfer rate from one 
quarter of the pipe may be written as the surn 

“ 9 P + 'll ■ »T> p . p P p + ff* L ) ( 33 ^ 

where the quarter-pipe effective perimeter is 

P p = ,R p /2 - b f S (34) 

Equation (33) is the basis for the optimization calculations in 
the following section. 

In this analysis the fin/pipe radiative interaction has been neg- 
lected for simplicity. This effect on overall heat rejection may be 
estimated from previous results (Sparrow and Cess, 1978), and a 
correction factor of 0.9 is used in Eq. (33); it is unlikely to have 
significant effect on the fin optimization. 

OPTIMIZATION PROCEDURE 

The optimization procedure used to determine the results given in 
this paper is presented here. The procedure was incorporated into 
a computer code, and several heat-pipe/fm combinations were 
analyzed. 

Determination of Fin Thickness and Length 


bj, respectively, and where the pipe is connected to the fin through a 
brazing, for example, with contact conductance, h c . Thus, the 
effective thermal resistance between the wick surface at vapor 
temperature (T = T v ) and the fin root can be written as 


1 - ifi ^ 

h e S k h c k p k l, 


(27) 


where S is a shape factor introduced to account for circumferential 
conduction in the pipe and liquid materials. The precise determi- 
nation of S requires a fully two-dimensional analysis of the com- 
posite region, which is beyond the scope of the present paper. It is 
noted, however, that the circumferentially constant temperature of 
the wick surface effectively limits the two-dimensional effect for heat 
flow into the fin, and we take S = 1.25 to 1.50. This contrasts with 
shape factors for convection pipes where the bulk fluid temperature 
is specified instead of a surface temperature (Lund, 1986, 1989). 


Because of this condition, there will be considerable direct heat 
rejection from the surface of the pipe at temperature T #p , 
effectively limiting the optimum fin length. The heat flux through 
the pipe wall is 



(28) 


Or, in terms of a pipe direct heat rejection efficiency 
«P“ , 'p e P <7T v 

where, from Eq. (28), the pipe efficiency is determined from 


7p - (1 " tip) 4 - 0 


(29) 

(30) 


and where the pipe radiation/resistance parameter is 


In order to determine the optimum fin geometry, the following 
in-puts are required: fin efficiency, surface Biot number, fin length, 
material thermal conductivities and densities, heat-pipe geometry 
and operating temperature. With this information, and specified S, 
Eq. (27) is used to calculate an effective thermal conductance h € 
between the inner surface of the heat-pipe wick and the fin root. 
Next, the effective end-conduction number 7 is calculated with 
Eq. (23). With 7 , fin efficiency and an estimated surface Biot 
number /?, the radiation conduction number N can be found from 
Fig. 4. Finally, the fin thickness b f is calculated using Eq. (16). 
This procedure is repeated for several different fin lengths until an 
optimum length is found. With the fin geometry fixed, a surface 
Biot number ft is calculated and compared with the initial esti- 
mate. If there is a significant difference between the two, the 
calculation is repeated with the new ft. However, in the cases 
considered, ft remained a very small number, on the order of 10 ' 4 , 
and iterations were unnecessary. 

Criterion Function 


The objective of the optimization procedure is to minimize the 
overall mass per thermal power dissipated (C, kg/kW), for the heat 
pipe and fin configuration. This gives a fin-tube arrangement which 
gives the maximum heat transfer per overall mass. One quarter of 
the mass per unit length of a heat-pipe and half of the fin is 

M = p[bfL + p p A p 4- />, A, (35) 

where pp p p , and P| are the densities of the fin material, pipe 
material and heat-pipe working fluid, respectively. A p and Aj are 
one-fourth the cross-sectional area of both heat-pipe material and 
heat-pipe working fluid, respectively. Half the total fin thickness is 
bp The total heat transfer from the heat pipe and fin (q t ) is calcu- 
lated using Eq. (33). The criterion function C, which represents 
the ratio of the mass of the heat-pipe and fin over the heat trans- 
fered, can then be calculated as 


5 



C = M/q t 


(36) 


.40 


The criterion function is calculated for a range of fin lengths 
and thicknesses generated from the calculation procedure discussed 
here. The fin length which yields the lowest C is chosen as 
optimum. Subsequently, conditions at this minimum, or optimum, 
are designated by a superscript asterisk; for example, the minimum 
value of the criterion is denoted C . 

Minimum-Mass Condition 


The procedure outlined in the preceding section to determine fin 
thickness and criterion function was incorporated into a computer 
code and used to model several different heat-pipe radiators. An 
aluminum heat pipe with aluminum fins and ammonia as a working 
fluid was modeled. Property data for this analysis are summarized 
in Table 1. Results from the analysis are shown in Fig. 6. The 
manner in which the criterion function varies with fin length for 
several different fin efficiency values is shown. It is interesting that 
the highest fin efficiency, rj f ~ 0.9, does not always yield the lowest 
criterion function C. The minimum value of C was calculated to be 
1.26 kg/kW. ThiB value corresponds to a fin efficiency of 0.5 and a 
fin length of 8.5 cm, which compares well with previous results with 
specified root temperature (Sparrow and Cess, 1978; Bartas and 
Sellers, 1960). There are a range of fin lengths and fin efficiency 
values which yield C values very close to the minimum of 1.26. 
Several combinations of parameters between fin lengths of approxi- 
mately 3 to 13 cm yield C values close to the minimum. This shows 
that in this particular case the radiator designer can use some 
discretion in choosing the optimum fin length, based on consid- 
erations other than thermal performance such as structural integ- 
rity, stowage volume, and survivability, etc. 


TABLE I. - INPUT PARAMETERS FOR 400 K 


RADIATOR MODEL 



Fin materia] 



Gr/AI 

Be/Al 

C-C/Al 

Working fluid 

nh 3 

nh 3 

nh 3 

Biot number Thermal 

8.1 xlO' 7 

3.4 xlO' 7 

6.3 xlO' 7 

conductivity, W/mK 




Fin length 

Varied 

180 

Varied 

Fin transverse 

77 

180 

10 

Pipe wall 

77 

180 

10 

Wick 

40 

40 

40 

Contact thermal 

Conductance between 


lx 10® 

1x10® 

pipe and fin, W/mK 

►— 1 1 

X 

o 

Material density, g/cm 3 




Fin 

2.702 

1.850 

1.69 

Pipe 

2.702 

1.850 

1.69 

Wick material 

2.702 

1.850 

2.80 

Heat pipe fluid 

0.37 

0.37 

0.37 

Wick porosity 

0.7 

0.7 

0.7 

Heat pipe o.d. cm 

1.5 

1.5 

1.5 

Shape factor, pipe/fin 

1.0 

1.0 

1.0 

Pipe wall thickness, mm 

1.0 

1.0 

1.0 

Wick thickness, mm 

1.0 

1.0 

1.0 

Emissivity 




Fin 

0.9 

0.9 

0.9 

Pipe 

0.9 

0.9 

0.9 

Vapor temperature 

400 K 

400 K 

400 K 

Sink temperature 

225 K 

225 K j 

225 K 

Radiation correction factor 

0.9 

0.9 ; 

0.9 



10 

Fin length, L, m 


Figure 6. — Criterion function versus fin length for several 
fin efficiencies (p ** 0; aluminum with ammonia heat-pipe 
fluid, geometry per table I). 


RESULTS 


Several heat-pipe and fin configurations made from different 
material combinations are modeled here. Metal matrix composites 
with graphite fibers, carbon-carbon, and beryllium are the fin and 
tube materials chosen for this analysis. The analysis is performed at 
two different vapor temperatures: 400 and 800 K. Two common 
heat-pipe working fluids are used; ammonia at 400 K and potassium 
at 800 K. Heat-pipe containment envelop materials are chosen 
which are compatible with the respective working fluid. The 
geometry of the heat-pipe is kept constant throughout the analysis. 
The results of this analysis yield an interesting comparison of some 
common space radiator materials. 

Graphite/aiuminum, carbon/carbon, and beryllium fin and tube 
configurations are designed here for operation at 400 K. Ammonia 
is used as the heat-pipe working fluid and aluminum as the wick and 
liner material for all three cases. Table I lists the input para- 
meters for this analysis. The results of the analysis are shown in 
Fig. 7. Two fin efficiencies ratios are shown for each case, = 0.5 
and = 0.9. These two values will typically bound the range of 
choices the radiator designer will have as shown in the optimization 
procedure section. Figure 7 shows a plot of the criterion function 
C versus axial thermal conductivity of the fin material. Axial ther- 
mal conductivity is varied because of the large variety of carbon 



Figure 7. — Criterion function as a function of fin axial thermal 
conductivity for three materials at two fin efficiencies 
ratios operating at a temperature of 400 K. 
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fibers available for use in the carbon or inetal matrix composites. 

The thermal conductivity of these fibers can vary from under 
10 W/mK to over 500 W/mK, However, the axial thermal conduc- 
tivity of beryllium is fixed. Therefore, the results for this case 
are represented by two points, as shown in Fig. 7. 

The results show that for a given axial thermal conductivity, 
carbon/carbon is superior to both beryllium and graphite/ 
aluminum. With P100 graphite fibers the axial thermal conduc- 
tivity of graphite/aluminum and carbon/carbon is calculated to be 
306 and 311 W/mK, respectively (Materials Science Corporation, 
1986). Specific thermal conductivities (axial thermal conductivity 
over density) of graphite/aluminum with P-100 fibers, carbon/ 
carbon with P-100 fibers, and beryllium are 121, 184, and 
97 (W/mK)/(g/cm 3 ), respectively. The specific thermal conduc- 
tivity of graphite aluminum with P-100 fibers is significantly greater 
than beryllium, however the beryllium yields a radiator design with 
a lower criterion function. This shows that materials with the 
highest specific thermal conductivity are not always the best 
material for a radiator application. This contrasts with earlier 
assumptions (McDanels et al., 1991). Low density appears to be 
more important than high thermal conductivity for this design. 

Also, the transverse thermal conductivity of the composite does not 
appear to significantly affect the criterion function. The transverse 
thermal conductivity of the carbon/carbon is significantly lower 
thau graphite/aluminum and beryllium as shown in Table l. How- 
ever, with both a low density and high axial thermal conductivity, 
the carbon/carbon is still shown to be superior to the other 
materials. 

A similar analysis was performed with graphite/copper, carbon/ 
carbon, and beryllium operating at 800 K. Potassium is the heat- 
pipe working fluid, and the wick and liner material is titanium for 
all three cases. The input parameters for this analysis are listed in 
Table II. The results of this analysis are shown in Fig. 8. Again, 


TABLE II. - INPUT PARAMETERS FOR 800 K 
RADIATOR MODEL 



Fin material 



Cr/Cu 

Be/Ti 

C-C/Ti 

Working fluid 

K 

K 

K 

Biot number 

2.4 xl0‘ 6 

4.8xl0' 6 

8.7 xl0‘ 5 

Thermal conductivity, W/mK 




Fin length 

Varied 

109 

Varied 

Fin transverse 

241 

109 

10 

Pipe wall 

241 

109 

10 

Wick 

7 

7 

7 

Contact thermal 

Conductance between 




pipe and fin, W/mK 

X 

o 

X 

H- > 

o 

lxlO 9 

Material density, g/cm 2 




Fin 

4.88 

1.850 

1.69 

Pipe 

4.88 

1.850 

1.69 

Wick material 

4.51 

4.51 

4.51 

Heat pipe fluid 

0.70 

0.70 

0.70 

Wick porosity 

0.7 

0.7 

0.7 

Heat pipe o.d., cm 

1.5 

1.5 

1.5 

Shape factor, pipe/fin 

1.0 

1.0 

1.0 

Pipe wall thickness, mm 

1.0 

1.0 

1.0 

Wick thickness, mm 

1.0 

1.0 

1.0 

Emissivity 




Fin 

0.9 

0.9 

0.9 

Pipe 

0.9 

0.9 

0.9 

Vapor temperature 

800 K 

800 K 

800 K 

Sink temperature 

225 K 

225 K 

225 K 

Radiation correction factor 

0.9 

0.9 

0.9 



Figure 8. — Optimum criterion function versus fin axial thermal 
conductivity for three materials at two fin efficiencies 
operating at a temperature of 800 K. 


the carbon/carbon design yields the lowest criterion function. The 
two beryllium cases are very close to the carbon/carbon case for the 
same axial thermal conductivity. The graphite/copper design yields 
a criterion function much higher than the beryllium case even 
though the specific thermal conductivity of graph! le/copper is 
greater than beryllium: 96 and 59 (W/mK)/(g/cm 3 ), respectively. 
This evidently results from the rather weak K x dependence of C , 
and the density proportional behavior of C in Eqs. (34) arid (35). 

It has been shown that the fin transverse conduction effect is 
negligible, except near the fin/pipe interface where two-dimensional 
conduction plays a role. This is incorporated in the shape factor S, 
where a larger value of S represents a larger two-dimensional effect 
in the fin root and heat-pjpe wall, and S = 1 represents one- 
dimensional radial conduction only. 

The effect of S on C is shown with an analysis of heat pipes 
and fins made of beryllium with potassium as the heat-pipe working 
fluid. The heat-pipe vapor temperature is 800 K. Figure 9 shows 
C is reduced with an increase in the two-dimensional (or circum- 
ferential) conduction effect. This is expected as a greater “group- 
ing” of energy at the root. Although beneficial, a shape factor in 
excess of i.5 is not expected for heat-pipe geometries as in Fig. 5, 
without special provision; exact evaluation of S requires a detailed, 
two-dimensional analysis of the interface region. Also evident in 
Fig. 9 is the^effect of fin efficiency on C , with 0.5 the value for 
minimum C . 



Figure 9. — Effect of shape factor on optimized performance 
index for beryllium/potassium heat-pipe/fin radiator panel. 
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CONCLUSION 

The analysis presented here shows that using specific thermal 
conductivity is not a good figure of merit for radiator materials. 

The results of the analysis showed that the density of a radiator 
material is much more critical than the axial thermal conductivity 
of that material in order to obtain maximum heat transferred with 
minimum mass. Also, the transverse thermal conductivity of a com- 
posite fin material is shown to have only a minor effect on the 
criterion functiou for the radiator designs analyzed. 
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